Properties of C in the ab initio nuclear shell-model 
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?— ( ! Abstract 
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\Q • We obtain properties of 12 C in the ab initio no-core nuclear shell-model. 

The effective Hamiltonians are derived microscopically from the realistic CD- 
. Bonn and the Argonne V8' nucleon-nucleon (NN) potentials as a function of 

the finite harmonic oscillator basis space. Binding energies, excitation spectra 

00 

\f~^ \ and electromagnetic properties are presented for model spaces up to 5hQ. The 

favorable comparison with available data is a consequence of the underlying 
■ NN interaction rather than a phenomenological fit. 

o : 
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+f ' While various methods have been developed to solve the three- and four-nucleon systems 

with realistic interactions few approaches are suitable for heavier nuclei at this time. 

Apart from the coupled cluster method |j] applied to closed-shell and near-closed shell nuclei, 
the Green's function Monte Carlo method is the only approach for which exact solutions of 

.£h ! systems with A < 8, interacting by realistic potentials, have been obtained ||. 

For more complex nuclei, treated as systems of nucleons interacting by realistic NN 
interactions, we apply the no-core shell-model approach [Brfl . To date, this ab initio method 



has been successfully applied to solve the three-nucleon as well as the four-nucleon bound- 
state problem ||||. Here, we address a vastly more complex system, 12 C, and present first 
results for an illustrative set of observables with two realistic NN interactions. 

There are several pressing reasons to investigate 12 C in a way that preserves as much 
predictive power as possible. The 12 C nucleus plays an important role |1C| in neutrino 



studies using liquid scintillator detectors. Also, there has been considerable interest recently 
in parity-violating electron scattering from (J 71 ", T) = (0 + , 0) targets, like 12 C, to measure 



the strangeness content of the nucleon |TT],|T2[. For these and many other reasons, there have 
been multi-frf2 shell model studies of 12 C in the past p!3|qT5|1. However, unlike our approach, 
phenomenological effective interactions were used. 

We start from the two-body Hamiltonian for the A-nucleon system, which depends on the 
intrinsic coordinates alone, Ha = T re i + V, where T re i is the relative kinetic energy operator 
and V is the sum of two-body nuclear and Coulomb interactions, V = Vn + Vc- There 
is no phenomenological one-body term. We neglect many-body interactions at present. To 
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facilitate our work, we add an A-nucleon Harmonic Oscillator (HO) Hamiltonian acting solely 
on the center-of-mass (CM), H CM = T CM + U C m, where U C m = ^Am^l 2 R 2 , R = \Ylt=i^i 
and m is the nucleon mass. The effect of this HO CM Hamiltonian will be subtracted in 
the final many-body calculation. The Hamiltonian, with a pseudo-dependence on Q, can be 
cast into the form 
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Since we solve the many-body problem in a finite HO model space, the realistic nuclear 
interaction in (|l|) will yield pathological results unless we derive a model-space dependent 
effective Hamiltonian. For this purpose, we adopt approaches presented by Lee and Suzuki 
1J ; 



Da Providencia and Shakin |17j, and Suzuki and Okamoto |L8|, which yield an Hermi- 
tian effective Hamiltonian. 

a unitary transformation of the Hamil 



According to Da Providencia and Shakin [17 



tonian H^, which is able to accommodate the short-range two-body correlations, can be 
introduced by choosing a two-body, in our case translationally invariant, antihermitian op- 
erator s = Y,£cj=i sucn 
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The transformed Hamiltonian can be expanded in terms of up to A-body clusters TC 
US 1 ' + + + . . . , where the one-body and two-body pieces are given as liS 1 ' 



Ef=i hi = E 



l<j: 



-i Va 



with 



V 12 = e- 5l2 (/ix + h 2 + V 12 )e Sl2 - (/i! + h 2 



(3) 



The full space is divided into a model or P-space, and a Q-space, using the projectors P and 
Q with P + Q — 1. It is then possible to determine the transformation operator S12 from 
the decoupling condition 



Q 2 e- 5l2 (/i! + h 2 + V 12 )e^P 2 = 



(4) 



The two-nucleon-state projectors {P 2 ,Q 2 ) follow from the definitions of the A-nucleon pro- 
jectors P, Q. This approach has a solution |18| , S\ 2 = arctanh(cj 
erator to satisfying u = Q 2 uoP 2 . This is the same operator, 
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with the op- 



which we previously em- 
ployed It can be directly obtained from the eigensolutions \k) of hi + h 2 + V\ 2 as 
(aQ\u\ap) = J2keK.( a Q\k) (k\oip), where we denote by tilde the inverted matrix of (ap\k). In 
the above relation, \ap) and |«q) are the 2-particle model-space and Q-space basis states, 
respectively, and JC denotes a set of dp eigenstates, whose properties are reproduced in the 
model space, with dp equal to the dimension of the model space. 

The resulting two-body effective interaction V\ 2 depends on A, on the HO frequency Q 
and on iV max , the maximum many-body HO excitation energy (above the lowest configu- 
ration) defining the P-space. It follows that TC^ + Ji^ 1 — Hqm is translationally invariant 
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and that Vu — > Vu for N mBX — > oo. A significant consequence of preserving translational 
invariance is the factorization of our wave function into a product of a CM |M1 component 
times an internal component. Hence, it is straightforward to correct exactly any observable 
for the CM effects. This feature distinguishes our approach from most phenomenological 
shell model studies that involve multiple HO shells. 

The most significant approximation used in the present application is the neglect of higher 
than two-body clusters in the unitary transformed Hamiltonian expansion. Our method 
is not a variational approach so the neglected clusters can contribute either positively or 
negatively to the binding energy. Indeed, we find that the character of the convergence 
depends on the choice of Q |||8||| • The method can be readily generalized in order to 
include, e.g., three-body clusters, and to demand the model-space decoupling on the three- 
body cluster level. Such a generalization leads to the derivation of the three-body effective 
interaction, which has been successfully applied in our calculations for the A = 4 system 
H0]. We learned from these studies that the contribution of higher order clusters to the 
ground-state energy are about 10% in similar model spaces that we employ here, when an 
optimal HO frequency is chosen, as it is the case in the present application. 

To solve for the properties of 12 C, we employ the m-scheme Many-Fermion Dynamics code 
T9| . Due to the fast growing matrix dimensions, reaching 6 488 004 at the iV max = 5 model 



space, we are restricted to iV max = 0, 2, 4 for the positive-parity states and iV max = 1, 3, 5 for 
the negative-parity states. Here, we utilize hfl = 15 MeV which lies in the range where the 
largest model space results are least sensitive to HQ. Full details will be reported elsewhere. 

We present results for the CD-Bonn |20j and the Argonne V8' Q NN potentials. Our 
positive-parity state results are presented in Table [I] and in Fig. [1], and the negative-parity 
state results are in Table [H] and Fig. |2|. While the energy of the lowest eigenstate of each 
parity increases with increasing model space, the relative level spacings are less dependent 
on model space size. As a gauge of trends with increasing model space size, consider the rms 
changes in excitation energies of the first 7 excited states of each parity in the CD-Bonn case. 
For positive parity states, the rms changes are 1.31 (0.22) MeV in going from to 2 (2 to 
A)hQ. For negative parity states, the rms changes are 0.87 (0.20) MeV in going from 1 to 3 
(3 to 5)M1 The difference between the iV max = 2(3) and 4(5) results is significantly smaller 
than that between the iV max = 0(1) and 2(3) results which is similar to the convergence 
trends we saw in lighter systems Our obtained binding energy of about 88 MeV in 

the ihfl space is expected to decrease with a further model space enlargement. We estimate, 
however, that our result should be within 10% of the exact solution for the two-body NN 
potential used. In order to reach the experimental binding energy, likely a true three-body 
NN interaction is neccessary M. 

In general, we obtain a reasonable agreement of the states dominated by Ohfl and lhfl 
configurations with experimental levels. We also observe a general trend of improvement 
with increasing model space size, i.e., the ordering of the T = 1 states. We obtain a resonable 
set of excitation energies for the T = 1 states relative to the lowest T = state of each 
parity. In addition, our lowest + T = 2 state lies between 27 and 29 MeV, depending on the 
NN potential and the model space, in good agreement with the experimental + 2 state at 
27.595 MeV. We note that the favorable comparison with available data is a consequence of 
the underlying NN interaction rather than a phenomenological fit. Our ground-state wave 
function in the 4hQ calculation contains about 61% of the OhQ component. The occupancy 
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of the 0p3/2 level is about 5.74 nucleons, while the occupancy of the Opl/2 level is about 
1.90 nucleons. From Tables | and [H], it is clear that the excitation energies of the negative- 
parity states relative to the positive-parity states decrease rapidly with the model-space 
enlargement. Still, even in our largest spaces the 3~0 state is more than 5 MeV too high 
compared to the experiment. 

In order to achieve a more realistic excitation energy a still larger HO expansion is needed 
especially for states with significant cluster structure. The two- and higher-/if2 dominated 
states, such as the 7.65 MeV + state that is known to be a three-alpha cluster resonance 
|?3fl , are not seen in the low-lying part of our calculated spectra. In general, the convergence 
rate of the 2hfl dominated states is quite different than that of the ground state as we 
observed in 4 He calculations performed in the present formalism [HJ§. Also, an optimal HO 
frequency for the convergence of the ground state will differ from the optimal frequency for 
the 2hQ states. We investigated the position of the lowest 2hQ dominated states and the 
giant-quadrupole resonance (GQR) E2 distribution. Our lowest 2hQ + state lies at about 
40 MeV excitation energy and the GQR E2 strength is fragmented between 43 to 50 MeV 
in the 2hfl calculation. In the 4M1 model space the excitation energy of the lowest 2hfl + 
state drops by 5 MeV to about 35 MeV and similarly the GQR strength position is lowered 
to 37-47 MeV. We note that the experimental GQR strength is observed in the range 18-28 
MeV [El. 



There is little difference between the results from the two NN interactions, although 
the overall agreement with experiment is slightly better for the CD-Bonn NN potential, in 
particular for the T = 1 states. 

Our radius and E2 results, based on the bare radius operator and bare nucleon charges, 
are smaller than the experimental values. The underestimation of the rms radius, the 
quadrupole moment and the E2 transitions is linked with the overestimation of the po- 
sition of the GQR strength and suggests that even in the N max = 4 space we still miss 
significant clustering effects. We also observe a strong model space dependence of the Ml 
transitions, 1 + 1 — > + 0. Clearly, there is still a need for effective operators, which are cal- 
culable within our theoretical framework. In general, to compute a two-body correction to 
a one-body operator in our formalism is more involved than the evaluation of the effective 
interaction. But, it is easy to study the lowest order renormalization for a two-body operator 
depending on the relative position of two nucleons as , e.g., the point-nucleon rms radius 
operator. Then, O e s ~ J2t<j=i e ~ Sij Oije Sij . We computed this term for the point-proton 
rms radius operator and found that the renormalization leads to an increase of the radius 
and that the size of this increase drops as the model space size increases. The r p results 
presented in Table | that were obtained without renormalization should be increased due to 
the renormalization by about 0.06, 0.02 and 0.01 fm for the iV max = 0, 2 and 4 model spaces, 
respectively. This does not imply that the renormalization of other operators, e.g., the E2 
operator, cannot be substantially higher. Similarly, as observed in our 3 H calculations j9], 
we anticipate that, in contrast with the energies, the higher-order corrections will be more 
significant and the overall convergence slower for other observables. 

We present these results as a useful description of the and lfrfi-dominated states of 
12 C. Our wave functions along with the one-body and two-body densities may also be used 
to predict cross sections for neutrino and muon reactions with 12 C. Such cross sections will 
be the subject of future investigations. The trends are encouraging and we will carry out 
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larger model space investigations in order to achieve greater convergence. 
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FIGURES 

FIG. 1. Experimental and theoretical positive-parity excitation spectra of 12 G Results obtained 
in AfiVL, 2%Q and O/jfi model spaces are compared. The effective interaction was derived from the 
CD-Bonn NN potential in a HO basis with h£l = 15 MeV. The experimental values are from Ref. 
& 

FIG. 2. Experimental and theoretical negative-parity spectra of 12 C. Results obtained in 5hft, 
3hQ and l?lO model spaces are compared. Other factors are the same as in Fig. [|. 
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TABLES 





12 C 




CD-Bonn 






AV8' 




model space 


- 


Ann 


2ftfi 


onn 


Ann 


2nn 


onn 


\E gs \ [MeV] 


92.162 


88.518 


92.353 


104.947 


87.675 


92.195 


104.753 


r p [fm] 


2.35(2) 


2.199 


2.228 


2.376 


2.202 


2.228 


2.376 


Q 2+ [e fm 2 ] 


+6(3) 


4.533 


4.430 


4.253 


4.536 


4.427 


4.250 


E x (0 + 0) [MeV] 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


£ x (2+0) [MeV] 


4.439 


3.697 


3.837 


3.734 


3.584 


3.766 


3.699 


eJi+0) [MeV] 


12.710 


14.141 


14.525 


13.866 


14.044 


14.549 


13.935 


£ x (4+0) [MeV] 


14.083 


13.355 


13.636 


12.406 


12.848 


13.255 


12.192 


£ x (l+l) [MeV] 


15.110 


16.165 


16.291 


15.290 


16.295 


16.515 


15.488 


£ x (2+l) [MeV] 


16.106 


17.717 


17.945 


15.970 


17.945 


17.823 


15.920 


£ x (0+l) [MeV] 


17.760 


16.618 


16.493 


14.698 


16.205 


16.208 


14.574 


B(E2;2+0 -> 0+0) 


7.59(42) 


4.625 


4.412 


4.092 


4.612 


4.397 


4.091 


B(Ml;l+0 -» 0+0) 


0.0145(21) 


0.0042 


0.0032 


0.0013 


0.0026 


0.0019 


0.0008 


B(Ml;l+0 -» 2+0) 


0.0081(14) 


0.0017 


0.0013 


0.0008 


0.0013 


0.0012 


0.0008 


B(M1;1+1 -» 0+0) 


0.951(20) 


0.355 


0.280 


0.158 


0.316 


0.252 


0.147 


B(M1;1+1 -» 2+0) 


0.068(9) 


0.0002 


0.0028 


0.0115 


0.0023 


0.0078 


0.0167 


B(E2;2+1 -> 0+0) 


0.65(13) 


0.283 


0.015 


0.0018 


0.104 


0.000 


0.002 



TABLE I. Experimental and calculated binding energies, ground-state point-proton rms radii, 
the 2^-state quadrupole moments, as well as E2 transitions, in e 2 fm 4 , and Ml transitions, in 
of 12 C. Results obtained in different model spaces, i.e., iV max = 4,2,0, and using effective 
interactions derived from the CD-Bonn and the Argonne V8' NN potentials are compared. A HO 
frequency nn 



15 MeV was employed. The experimental values are from Refs. pl,B2 
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12 C 




CD-Bonn 






AV8' 




model space 


- 




3M2 


inn 




3nn 


mn 


|^(3-0)| [MeV] 


82.521 


72.952 


75.331 


83.390 


72.300 


75.360 


83.459 


r p [fm] 




2.309 


2.316 


2.425 


2.310 


2.315 


2.425 


Q 3 - [e fm 2 ] 




-7.942 


-7.596 


-6.936 


-7.920 


-7.575 


-6.933 


E{3~0) - E gs [MeV] 


9.641 


15.566 


17.022 


21.557 


15.375 


16.835 


21.294 


£ x (3"0) [MeV] 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


eJi-0) [MeV] 


1.203 


2.093 


2.256 


1.561 


2.112 


2.274 


1.552 


E x (2-0) [MeV] 


2.187 


3.722 


4.051 


3.582 


3.722 


4.057 


3.567 


E x (4-0) [MeV] 


3.711 


4.866 


5.084 


4.768 


4.741 


4.993 


4.710 


E x (0-0) [MeV] 




7.148 


7.062 


5.712 


7.148 


7.156 


5.777 


^(2-1) [MeV] 


6.929 


7.671 


7.783 


7.340 


7.949 


8.237 


7.574 


EJ3-0) [MeV] 




7.877 


8.151 


6.886 


7.651 


7.983 


6.745 


£ x (l - l) [MeV] 


7.589 


8.048 


7.951 


7.042 


8.117 


8.096 


7.184 



TABLE II. Experimental and calculated negative-parity state energies, the 3~0-state 
point-proton rms radii, and quadrupole moments are shown. Results obtained in different model 
spaces, i.e., N max = 5,3,1, and using effective interactions derived from the CD-Bonn and the 
Argonne V8' NN potentials are compared. The calculated excitation energy of 3~0 is obtained 
by comparing its energy in the NhQ space with the ground state in the (JV — l)7iQ space. A HO 
frequency h£l = 15 MeV was employed. The experimental values are taken from Ref. [21]. 
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10 
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12 



c 



^Q=15 MeV 



Exp 




o + o 

2 + 1 
1 + 1 
+ 1 
2 + 

1 + 
4 + 



2 + 



+ 



4*Q 2^Q 



6 



CD 



8 



2" 1 



^ 4 

fcq 4- o 



2" 

r o 



12 



Q 1iQ=15 MeV 



5 2- l 

r l • . " / 




i- l 

3" o 

o- o 
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2" 
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L 3" 3" 

Exp 5^Q 3^Q VhCl 



